Double zigzag spin chain in strong magnetic field close to saturation 
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We study the ground state phase diagram of a frustrated spin tube in a strong external magnetic 
field. This model can be viewed as two coupled zigzag spin chains, or as a two-leg spin ladder with 
frustrating next-nearest-neighbor couplings along the legs, and its study is motivated by the physics 
of such materials as Sulfolane-Cu2Cl4 and BiCu2P06. In magnetic fields right below the saturation, 
the system can be effectively represented as a dilute gas of two species of bosonic quasiparticles that 
correspond to magnons with inequivalent incommensurate momenta at two degenerate minima of 
the magnon dispersion. Using the method previously proposed and tested for frustrated spin chains, 
we calculate effective interactions in this two-component Bose gas. On this basis, we establish the 
phase diagram of nearly-saturated frustrated spin tube, which is shown to include the two-component 
Luttinger liquid, two types of vector chiral phases, and phases whose physics is determined by the 
presence of bound magnons. We study the phase diagram of the model numerically by means of 
the density matrix renormalization group technique, and find a good agreement with our analytical 
predictions. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Mg, 67.85.Hj 



I. INTRODUCTION 



Frustrated spin systems, especially in low dimensions, 
display a rich variety of unconventionally ordered ground 
statesjii^ Strong external magnetic field, competing with 
the exchange interaction, can serve as a control param- 
eter that drives the corresponding quantum phase tran- 
sitions. The ground state of a frustrated quantum spin 
system is considerably simplified in a sufficiently strong 
external field that eventually leads to a fully polarized 
state above some critical field value Hs (strictly speaking, 
the latter is true only for an axially-symmetric case, but 
we assume that deviations from axial symmetry are neg- 
ligibly small). In fields just slightly below Hg, one may 
view the system as a dilute gas of excitations (magnons) 
on top of the fully polarized statci^"^ At low density of 
magnons they can be approximately treated as bosonic 
quasiparticles. In the case of a strong frustration, the 
magnon dispersion has two or more degenerate minima 
at inequivalent incommensurate wave vectors, so one ar- 
rives at the picture of a multicomponent dilute Bose gas. 
In the one-dimensional case, infrared singularities, ap- 
pearing in the description of effective interactions in the 
magnon gas, require special treatment 

Depending on the ratio of interactions between the 
same or different sorts of particles, several types of the 
ground state can be favored. Particularly, in two- and 
three-dimensional systems, different kinds of helical or- 
der ("fan" and "umbrella") are realized)^'^ while in one 
dimension quantum fluctuations destroy long-range heli- 
cal order and may lead to the formation of several differ- 
ent states with competing types of unconventional short- 
and long-range orders. In one dimension, in the case of 



repulsion between magnons, the "umbrella" and "fan" 
phases get replaced by the vector chiral (VC) long-range 
order— ^— (which is equivalent to the local spin current) 
and by the two-component Tomonaga-Luttinger liquid 
(TLL2)^ respectively. On the other hand, attraction 
between quasiparticles can lead to the appearance of a 
short-range multipolar (spin nematic) order^iii^ or al- 
ternatively to metamagnetic jumpsJ^ 

Recently, the above approach, based on the mapping to 
the multicomponent Bose gas, has been successfully ap- 
plied to spin-^ zigzag spin chain^ii^ which is a paradig- 
matic model of a frustrated spin system. It has been 
shown that for zigzag chains close to saturation this ap- 
proach is able to capture the physics of phase transitions 
between the VC and TLL2 phases, and for 5' > 1 it can 
detect the boundary of the (metamagnetic) region where 
bound states of magnons are formed. In the present pa- 
per, we employ this method to study the strong-field part 
of the ground state phase diagram of the frustrated spin 
tube shown in Fig. [T] The spin tube, which will be the 
subject of our study, can be viewed as two coupled zigzag 
spin chains, or as a two-leg spin ladder with frustrating 
next-nearest-neighbor couplings along the legs, see Fig. 
[T] This model is described by the following Hamiltonian: 

L 

n=l a=l,2 

where Sn.a are spin-S* operators acting at the rt-th site of 
the a-th leg, Ji and J2 are the nearest-neighbor (NN) and 
next-nearest neighbor (NNN) exchange couplings along 
the legs, J is the rung exchange, and H is the external 
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FIG. 1. (Color online). The frustrated spin tube described 
by the Hamiltonian ([1} . The tube can be alternatively viewed 
as two zigzag chains coupled by the transversal interaction 
J_i, or as a spin ladder with next-nearest-neighbor exchange 
couplings along the legs (lower panel). 



magnetic field. The system may be alternatively viewed 
as four antiferromagnetic chains connected by rung and 
zigzag couplings. This model has been recently stud- 
ied at zero fieldj ^"'^^ It is believed to be relevant for 
the physics of such quasi-one-dimensional materials as 
Sulfolane-Cu2C;4 (Cu2Cl4H8C4S02),22 which exhibited 
unusual critical behavior in the field-induced transition 
from a helimagnetic to a non-magnetic phase, ^'^^^^ and 
BiCu2P06, which has attracted the attention of several 
research groups as being a realization of a frustrated lad- 
der system with incommensurate correlations. ^^-^ Apart 
from the possible relevance for the above materials, this 
model is fundamentally interesting since it presents the 
simplest example of two interacting zigzag chains. 

In this paper, we are interested in the frustrated case, 
so J2 is chosen to be positive while Ji may have any sign. 
It is convenient to use the quantity 



/3 = J1/J2 



(2) 



as the frustration parameter. We will be interested in 
the regime |/3| < 4, when the magnon dispersion develops 
two degenerate minima at inequivalent points ±Q in the 
momentum space (in what follows, we refer to this regime 
as "strong frustration" ) . 

For S* = i, the phase diagram of the above model 
in the absence of the magnetic field has been studied 
numericallySi and was found to contain the rung singlet 
and the columnar dimer phase. Earlier, a slightly dif- 
ferent version of the model including exchange coupling 
along diagonals of the ladder has been investigated^^ its 
phase diagram at zero field has been shown to contain 
the rung singlet phase, the Haldane phase, and two dif- 
ferent columnar dimerized phases. The magnetic phase 
diagram of both versions of the model is at present un- 
explored. 

We study the ground state of the strongly frustrated 
(|/3| < 4) spin-S* tube decribed by the Hamiltonian ([T]) in 
high magnetic fields in the immediate vicinity of satura- 
tion, for spin values S = 1 and S = ^. It is shown that 



the phase diagram contains the two-component Luttinger 
liquid, two types of vector chiral phases, and phases 
whose physics is determined by the presence of bound 
magnons. We compare our analytical predictions with 
the results of numerical simulations using the density ma- 
trix renormalization group ^^'^° (DMRG) technique. To 
that end, we compute the chirality correlation function 
and magnetization distribution as functions of /3 at sev- 
eral values of J_l , at fixed magnetization close to satura- 
tion. We demonstrate that the DMRG results are in a 
good agreement with our theoretical predictions. 

The structure of the paper is as follows: in Sect. [Tl] 
we describe the mapping of the spin tube problem to the 
dilute two-component lattice Bose gas and outline the 
main steps of computing effective interactions. Section 
mil discusses the specific predictions of the theory for spin 
tubes with S = \ and 1, while Sect. IIVI presents the 
results of numerical analysis and their comparison with 
analytical predictions. Finally, Sect. |V] contains a brief 
summary. 



II. EFFECTIVE TWO-COMPONENT BOSE 
GAS DESCRIPTION OF THE SPIN TUBE 

We intend to map the spin problem ([1]) to a dilute gas 
of interacting magnons, for values of the field H slightly 
lower than the saturation field Hs- For that purpose, it 
is convenient to use the Dyson-Maleev representation for 
the spin operators in ([l}: 



= V2Sh 

-s-bi 



S-nrn — \/2S'&J,„, ( 1 



bi 



2S 



(3) 



where bnm are bosonic operators acting at site (nm) of 
the lattice, and n = 1, . . . , L and m = 1, 2 denote the 
rung and leg numbers, respectively, see Fig.[T] To enforce 
the constraint &Jim&„m < SS*, one can add the infinite 
interaction term to the Hamiltonian, which reads: 



n^n + uY,:{blj>^j 



2S+1 



u 



(4) 



where : (...) : denotes normal ordering. At the level of 
two-body interactions (which are dominating because of 
the diluteness of the gas) this term should be taken into 
account only for S — ^. 

Assuming periodic boundary conditions, we pass to the 
momentum representation for bosonic operators. 



bnm = b^e 



where k — {kx,ky), with ky taking only values or tt, 
and L is the total number of rungs. Then one can cast 
the Hamiltonian ([T]) in the following form: 

n = E ^'^^1^-^ + 27 E Uk,k')bl^^bl_^b^b^,. (5) 



kk'q 
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Here the magnon dispersion Ef. is given by 

Ek=H- 2(Ji + J2)S + 2SJk^ + J^S{cosky - 1), (6) 

where we use the shorthand notation 

Jfc = Ji cos(fc) + J2 cos(2fc). (7) 

In the case of the strong frustration \/3\ < 4, which is 
of main interest for us, there are two inequivalent de- 
generate minima of that are reached at wave vectors 
k = {±Q,n) and k = (±(5,0) for positive and negative 
J±, respectively. The wave vector Q is incommensurate 
and is given by 

Q = arccos(-^/4). (8) 

The saturation field Hg can be found^ from the condition 
(min Ek)\fj^^^ = : 

Hs=2S{Ji + J2 + J±e{J^)-jQ}, (9) 

where 6{x) is the Heaviside function. The external field 
may be viewed as playing the role of the chemical po- 
tential ^ = Hg — H ioi magnons. In what follows, it is 
convenient to introduce instead of Ek the quantity 

ek = Ek + fi. (10) 
I 



For 5* > 1, the two-body interaction Vq{k, k') depends 
on the transferred momentum q as well as on the incom- 
ing momenta k, k': 

Vqik,k')=2Jg^-Jk^-Jk'^ (11) 

-I- J±{cosqy - ^{cosky + cosk'y)}. 

For spin i, one has to add the term ^ to the Hamil- 
tonian, simultaneously dropping the terms like hnmbnm 
involving double occupancy. As a result, for S = ^ the 
expression for the two-body interaction simplifies to 

Vq{k,k') = U + J^CQsqy + 2Jq^, [/ ^ +00. (12) 

The model ([5]) describes two magnon branches of differ- 
ent parity with respect to the permutation of the ladder 
legs. We denote the operators describing even and odd 
magnons by Cfc and Ofc, respectively: 

afe=&(fc,7r): Cfc=6(fco)- (13) 

In this notation, the Hamiltonian takes the form 



n = Y,[{el- m)4«. + (4 - m)4c,} + — ^ {l/-(fc, k')al^^al,_^a,a^, + ^-(fc, fc')4+,4,_^c,c,, 

k kk' q 

+ Vq^^{k,k')al^^al,_^c^c^,+V^'^{k,k')cl^/^,_^a^^^^ 



The magnon energies and interaction amplitudes above 
can be read off Eqs. (HI), ([10]), (HI]), For the energies, 
one has 

£^ = 25(Jg- Jq), e^ = e^ + 2Ji5 at > 0, (14) 
e'i^ = 2S{Jq~ Jq), £^ = e^ + 2|Ji|5 at < 0, 

thus the energies of a-branch lie below (above) those of 
the c- branch for > ( J_i < 0), respectively. When the 
magnetic field is decreased below iJ^, the ground state of 
the system can be viewed as a dilute gas of a-magnons 
for Jj^ > or of c-magnons for J± < Q. Therefore, the 
last term in (fH)) . describing the scattering of a c- magnon 
on an a-magnon, does not influence the structure of the 
ground state in the immediate vicinity of the saturation 
field: under the condition 

Hs - H ^ Ai = 2SJ± (15) 

(see Fig. [2]) there is simply no regime when densities 
of both sorts of magnons (a and c) are simultaneously 
nonzero. In what follows, we assume that the condi- 
tion (fT5)) is always satisfied, so we can safely ignore the 



presence of the last term in (fT4)) . However, the other 
amplitudes in (|14p. e.g., describing conversion of a pair 
of a-magnons to a pair of c-magnons, have to be kept, 
because they contribute to intermediate virtual states in 
multiple scattering processes. 

The other two-body interaction amplitudes for S" > 1 
are 

y-(fc, k') = y™(A:, k') i(jfc + jfcO, 

F;-(fc, k') ^J,-J^- i(Jfc + Jk'), (16) 

Vq^''{k,k') = J, + - ^{Jk + Jk'), 
and for S — ^ they have to be modified as 

V-c = V^'^ = \{U-J^) + J,. (17) 

We have thus mapped the initial spin problem onto 
a Id lattice gas of particles with a nontrivial double- 
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FIG. 2. (Color online). Schematic picture of the magnon 
dispersion in the vicinity of the saturation field. For the sake 
of definiteness, the case of antiferromagnetic J± is shown. The 
system is populated by a-magnons with the momenta close to 
±Q, provided that the chemical potential ^ — Hs — H satisfies 
the conditions (|15|) and (|19|) . 



minima dispersion. The renormalized two-body interac- 
tion in such a gas can be easily found in the dilute limit, 
i.e., /Lt — >■ 0. Since the self-energy vanishes at /x — 0, 
the full propagator coincides with the bare one^^ and 
thus the Bethe-Salpeter (BS) equation for the renormal- 
ized two-body interaction vertex r^'^(fc, fc'; E) (where E 
is the total energy of the incoming particles) takes the 
following formi^ 



rf (fc,fc';i?) = C'3(fc,fc') 



(18) 



L 



EE 



V^2^{k+p,k' -p)Tf{k,k'-E) 



'-k+p 



E 



where labels a, /?, 7 denote the magnon branch and can 
take the values "a" and "c". The above equation is 
schematically shown in terms of Feynman diagrams in 
Fig. ®. 

If the magnetic field is close enough to the saturation, 
so that the condition 



p = Hs-H<^A2 = 2J2S{\I3\/A - 1)2 



(19) 



is satisfied (see Fig. then the system is mainly pop- 
ulated by magnons (of a- or c-branch, depending on the 
sign of Jl) with momenta around the two dispersion 
minima at itQ, which at low energies can be interpreted 
as two different bosonic "flavors". For those low-energy 
modes, one can formulate the effective theory in the form 
of the Gross-Pitaevsky-type energy functional for a two- 
component Bose field: 



n 



GP 



/^-{E 



|V$ 



CT|2 



2m 



Ti2nin2 



^{ni + 712) I ■ 



(20) 



Here the Planck constant is set to unity, are the 

macroscopic bosonic fields that describe magnons with 
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FIG. 3. The ladder approximation to the Bethe-Salpeter 
equation for the renormalized two-body interaction vertex 
Tq°'{k,k' ; E). Solid lines denote bare propagators. The ap- 
proximation becomes exact at ^ — >■ 0. 



momenta k lying within the intervals |fc±(5| < A around 
the dispersion minima, A is the infrared cutoff, — 
l^'^P are the particle densities, and m is the effective 
mass, 



1 

TO 



5J2(16-/?2) 



(21) 



The point Fn — V12 corresponds to the enhanced SU{2) 
symmetry at the level of the effective low-energy the- 
ory. For < Fii, the ground state of the gas con- 
tains an equal density of the two particle species, and for 
> Fii just one of the two species is present in the 
ground state. In the spin problem the total number of 
each bosonic species is not separately fixed, in contrast to 
a typical setup for atomic mixtures. In a setup with fixed 
particle numbers, the ground state at < Fn is in the 
mixed phase, and F12 > Fn corresponds to phase sepa- 
ration. In the spin language, the separated phase maps 
to the vector chiral (VC) phase, while the mixed phase 
corresponds to the two-component Tomonaga-Luttinger 
liquid (TLL2)J^»^ 

Macroscopic effective couplings Fn, F12 of the Gross- 
Pitaevsky-type theory can be obtained from the E ^ 
limit of the corresponding vertex functions: 



rii = Fr(Q,Q;0)U^^^, 
ri2= [Fr(-Q,Q;0) 



"paa / 
^ 2Q[ 



-Q,Q;0) 



lA=A, 



(22) 



■'c" for positive and negative J_l, re- 
spectively, and the vertex function T'^°'{k,k'\E) is the 
solution of the BS equation (fT8)) with the infrared cut- 
off \p\ > A employed in the summation over internal 
transferred momenta p. The resulting expressions are 
viewed as functions of the running cutoff A in the spirit 
of the renormalization group (RG) approach, and the 
RG flow A is then interrupted at a certain scale 
A = A* = pm/2 that depends on the chemical poten- 
tial (or, in other words, on the magnon density). The 
above approach is well known for one-component Bose 
ga3^"~ and has been successfully applied to the multi- 
component case recently.— 

There is an alternative approach ; ^"'^^i'^" which, instead 
of the infrared cutoff in the momentum space, introduces 
an "off-shell" regularization: the two-body scattering 
amplitudes in the presence of a finite particle density are 
obtained by taking the "bare" expressions V^^ (k, k'; E) 
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at a finite negative energy E = — i?* = — 7r^/i/8: 

rn = rr(Q,0;-£^*)lA^o' (23) 
ri2 = [rr(-Q, Q; -e.) + r^§(-g, Q; 1^^^ , 

where a takes the same value as in Eq. ((22|) . One can 
show-'^^ that the off-sheh regularization yields the results 
that are equivalent to the cutoff regularization scheme. 
In this work, we have used the off-shell regularization 
because it is more convenient technically. 

Our model contains only short-range interactions, so 
the solution of (IT5| can be expressed in terms of a finite 
number of Fourier modes in the transferred momentum<^ 
In our case, from the structure of V^^{k,k') it is easy 
to see that each component of r^''(fc, fc'; E) can contain 
only five Fourier harmonics proportional to 1, cosg, sing, 
cos 2g, and sin 2q. The system of integral equations is 
thus reduced to a system of linear equations that can be 
solved analytically for any value of the spin S. 

For the purpose of finding only the Gross-Pitaevsky ef- 
fective couplings Fii, Ti2, the problem can be simplified 
even further. First of all, the system (|18l) describes four 
equations for the vertices, which split into two decoupled 
pairs: a pair of coupled equations for F"", F™, and an- 
other pair of coupled equations for F^^ and F"'^. For Jj^ > 
0, the lowest energy excitations are a-magnons, and thus, 
in order to find the effective couplings Fn, F12, we only 
need to solve the first pair of the BS equations for F°°, 
F°"; similarly, for Jj^ < we are interested only in the 
equations for F'^'^ and F°^. Second, it is easy to see that 
Fii can be found as the q — Q value of F^"((5, Q; —E^) 
which is an even function of the transferred momentum 
q, and F12 may be represented as the q — Q value of the 
function T'^1^y{-Q,Q;-E,) + r2Jtg{-Q,Q;-E,), which 
is also even in q. For that reason, one can rewrite the 
integral equations (fT8|) for the above two even functions, 
keeping only even Fourier harmonicsi^ This reduces the 
number of resulting linear equations to six and makes the 
problem amenable to analytical treatment. We refer the 
reader to the Appendix for further details. 

Solving the Bethe-Salpeter equation, one can show 
that the expansion of Fn, F12 in i?* has the following 
structure>iS 

Fii \4:E^J gii 

- TTT + — + 0{EI'^) + ... 24 

F12 \iE^ J gi2 

Note that for E„ ^ (i.e., ^ — ?> 0) the effective couplings 
Fii and F12 flow to the same value, which reflects the 
tendency of the RG flow to restore the SU (2) symmetry 
for the two-component Bose mixture.'^* 

Parameters 511, gi2, entering the second term in the ex- 
pansions ([M)) . under certain conditions, namely |5ii|m ^ 
1 and |(/i2|w ^ 1, can be identified with the effective 
bare coupling constants of the continuum two-component 
Bose gas with contact interactions (see Ref. ^id for de- 
tails). If the above conditions are broken, parameters 



gij cannot be interpreted as physical bare couplings, and 
only the renormalized interactions Fy retain their mean- 
ing as effective low-energy coupling constants. The only 
physical meaning of gij in such a case is that they are con- 
nected to the asymptotic phase shift of scattering states 
at small transferred momentaJ^ 

From Eqs. (IMl) one can see that transition points be- 
tween the TLL2 and VC phases, that are determined by 
the condition Fn = F12, correspond only to crossings of 
512 and gii. (See, for example. Fig. 2^: when gu goes 
through a pole changing sign from plus to minus infinity, 
then at the pole gi2 — gn changes sign from negative to 
positive, but this does not correspond to any phase tran- 
sition, because on both sides of the pole in its immediate 
vicinity F12 < Fii). 

Similarly, gn or 1712 becoming negative by going 
through a zero indicates the appearance of magnon 
bound states with zero total momentum, while a change 
of sign through a pole is not signaling any transition, 
but rather indicates a crossover into the so-called "super- 
Tonks" regime. It should be emphasized that within 
the present effective theory, which is essentially based on 
the two-body interaction, we cannot predict whether the 
formation of bound states stops at the level of bound 
magnon pairs, or continues with multiparticlc bound 
states. 



III. STRONG-FIELD PHASE DIAGRAM: 
ANALYTICAL RESULTS 

Let us turn our attention to the specific predictions of 
our theory for the frustrated spin tube model defined by 
the Hamiltonian ([IJ , at two spin values S = 1 and 5 = ^ . 
Details concerning the solution of the BS equations can 
be found in the Appendix. 

Fig. |4^-c illustrates the behavior of the "bare" effec- 
tive couplings 511, gi2 for the spin-1 tube, as functions of 
the frustration parameter j3, for several values of the in- 
terchain coupling J±. VC-TLL2 transitions are detected 
by crossings of 171 1 and gi2, and zeros of 51 1 signal the 
formation of magnon bound states with the total momen- 
tum k — ±2Q. The resulting phase diagram is shown in 
Fig. HJi. One can see that the region of small |/3| always 
corresponds to the chiral phase, similar to the case of a 
single frustrated spin-1 chain. ^"^'^^ For antiferromagnetic 
zigzag coupling /3 > 0, there is only one transition be- 
tween the vector chiral and the two-component Luttinger 
liquid phases, which is rather weakly dependent on the 
interchain (rung) coupling J±. Nonanalytic behavior of 
the phase boundary at = stems from the follow- 
ing: we assume that we work in the immediate vicinity 
of the saturation field (see conditions (|15l) and (ITOl) ). the 
ground state contains only a-magnons at Jj_ > and 
only c-magnons at J_l < 0. It is clear that the magnetic 
field range, where the assumption (|15p remains applica- 
ble, shrinks to zero as J± 0, which causes the above 
nonanalyticity. At Jj^ = the model corresponds to two 
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FIG. 4. (Color online), (a) "Bare" coupling constants gu, gi2 of the eflFective theory, obtained from the analytical solution 
of regularized Bethe-Salpeter equations, for 5* = 1 frustrated spin tube with J±/J2 = 2; (b), (c) the same for J_l/J2 — 0.02 
and J±/J2 ~ —2; (d) the predicted phase diagram in the vicinity of saturation; symbols show the transition points obtained 
from numerical simulations. Crossings of gn and gi2 correspond to transitions between the vector chiral and two-component 
Luttinger liquid phases, while the regions where gn becomes negative by going through a zero (not through a pole) indicate 
the appearance of magnon bound states. 



decoupled frustrated chains, so two magnon branches be- 
come degenerate, and the problem reduces to that for a 
single chain.— 

For ferromagnetic zigzag coupling /3 < 0, there is a 
large region with negative gn supporting bound magnon 
states. From the numerical analysis for a single S* = 1 
frustrated chain^ it is known that at least around J± = 
there is a metamagnetic jump in the magnetization 
curves in this region. Presence of such a jump indicates 
formation of "magnon drops" - bound states of a large 
number of magnons. Increasing antiferromagnetic rung 
interaction J± > leads to the opening of a finite TLL2 
phase window close to /3 = —4 (which is the boundary of 
a transition into a one-component Luttinger liquid state) . 

Fig. [S] shows two examples of the characteristic behav- 
ior of the bare coupling, along with the resulting phase 
diagram, for the S — ^ tube. The topology of the phase 
diagram is qualitatively the same as in the S — 1 case, 
but there are certain caveats which one should have in 
mind. One important difference concerns the region of 
bound magnon states. It is known^ that for a single 
isotropic •S' = ^ frustrated chain with /3 < the lowest 
energy of a two-magnon bound state is not reached at 



the total momentum k — ±2(5, as one could expect from 
the picture of two bound magnons of the same flavor, 
but instead the minimum of the bound state dispersion 
lies at fc = TT in a rather wide region of —2.67 < /3 < 0. 
For that reason, one may expect that the actual size of 
the region dominated by magnon bound states is larger 
than that of the region labeled "fc = 2Q bound states" 
in Fig. [21 Second, from the analysis of the single-chain 
S = ^ problem in Ref. [l^, it is known that predictions of 
the present theory for the VC-TLL2 transition at Jj_ = 
should deviate substantially from the numerical results. 
In Section llVI we will see, however, that agreement with 
the numerics is actually improved with the increase of 
the rung coupling strength | J_l | . 



IV. NUMERICAL ANALYSIS 

To verify our analytical predictions, we have studied 
the frustrated spin tube model ([T]) with S — 1 and 
S — ■T^ using the density matrix renormalization group^^ 
(DMRG) method (see Ref. [s^ for a detailed description 
of the DMRG technique). 
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FIG. 5. (Color online). "Bare" coupling constants gn, (;i2 
for S — \ frustrated spin tube with (a) J^j Ji = 2 and (b) 
J±/ J2 ~ —2; (c) the predicted phase diagram in the vicinity 
of saturation, symbols show the transition points obtained 
from numerical simulations. 



We study the vector chirality correlation functions to 
identify the phases that have long-range vector chiral or- 
der. The ground state of spin tubes with J± > near 
the saturation is populated with the magnons of the an- 
tisymmetric a-branch, so the relevant quantity in that 
case is the antisymmetric chirality 



I^A{n) = {Sn,l - Sn,2) X {Sn+1,1 - S'„+i, 



(25) 



Similarly, for J± < Q one has to look at the correlators 
of the symmetric chirality 



(26) 



In order to identify regions with metamagnetic be- 
havior, i.e., the regions where magnon attraction leads 
to the formation of a single bound state consisting of 
a macroscopic number of magnons ( "magnon drop" ) , 
we calculate the distribution of the rung magnetization 
Mn = {S^ 1 + 3^2) along the tube. This approach, how- 
ever, does not allow to detect phases where the formation 
of bound states stops at the level of a finite number of 
magnons. 

We use DMRG in its matrix product state 
formulation)^!^ which allows us to exploit the non- 
Abelian SU{2) symmetry, as well as the Abelian U{1). 
(While the magnetic field H breaks SU{2) symmetry, 
the fact that the Zeeman energy term commutes with 
the rest of the Hamiltonian makes it possible to take the 
influence of the magnetic field into account by calculat- 
ing the ground state of the model in a sector with the 
given total spin Stot-) The advantage of using the SU{2) 
symmetry lies in a considerable reduction of the number 
of states TO which is necessary to describe the system, 
because one essentially treats the multiplet of states of 
the same total spin as a single representative state. 

The use of SU{2) symmetry has a disadvantage as well: 
since the non- Abelian method allows to compute only re- 
duced matrix elements (in the sense of the Wigner-Eckart 
theorem), one can only compute rotationally invariant 
correlators such as {KA{n) ■ KA{n')), etc. This can be 
inconvenient if the contribution of the transversal com- 
ponents of chirality exhibits strong oscillations that act 
as a "noise" masking the long-range order in the longi- 
tudinal component, as it has been found in frustrated 
chains4^ We have found such strong oscillations for ks 
in spin tubes with J± < 0, while the correlators of ka 
in Jj^ > systems were essentially free from oscillations, 
calculating 

For that reason, we have used SU (2) symmetry in 
our calculations for antiferromagnetically coupled tubes 
{J± > 0), and resorted to the standard t7(l) calculations 
for the J± < case. Fig. [5] shows typical examples of 
chiral correlators for systems with ferro- and antiferro- 
magnetic sign of J±, calculated with or without the use 
of the SU (2) symmetry. 

We have studied spin-1 and spin-i spin tubes con- 
sisting of up to 256 spins, with open boundary condi- 
tions, keeping up to m = 600 states in most calculations. 
The total magnetization M of the system has been set 
at about 90% of the saturation value Ms (specifically, 
we have kept M/Ms = 115/128 for S" = 1 tubes and 
M/Ms = 29/32 for S" = i ones). The (squared) chiral 
order parameters k\, k| were extracted from the large- 
distance behavior of the corresponding correlation func- 
tions (the technicalities of this procedure are described 
in detail in Ref. 03). 

Fig. [7] shows the behavior of the chiral order param- 
eters along three constant- Ji cuts in the phase dia- 
gram of the spin-1 tube. For antiferromagnetic rungs 
(J_l/J2 = 2) one can clearly see two transitions around 
/? w 2.1 and /3 ss —3.40, which is consistent with the pre- 
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FIG. 6. (Color online). Typical DMRG results for the chiral- 
ity correlators of the L — 64 spin-1 tube (128 spins) at the 
total magnetization Stot ~ 116, inside the vector chiral phase: 
(a) rotationally invariant correlator of the antisymmetric chi- 
rality {KA{n) ■ ka{0)), calculated with the use of the SU{2) 
symmetry: there are no visible oscillations; (b) rotationally 
invariant correlator of the symmetric chirality {Ks{n) ■ ks{0)) 
exhibits strong oscillations which are due to the contribution 
of the transversal components as seen from (c) longitudinal 
correlators of the same quantity. 
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FIG. 7. (Color online). Chirality order parameters of the 
5=1 tube obtained from the large-distance behavior of the 
correlation functions. The results for L = 64 tube (128 spins) 
were calculated in the sector with the total spin Stot = 116, 
and for L = 128 we took Stot = 230. (a) J±/J2 = 0.5; (b) 
J±/J2 = 2.0, only L = 64 results are shown; (c) J±/J2 = 
—2.0, here the calculations have been done using the usual 
U{1) DMRG method for L = 128 system. The point at ^ = 
has not been obtained numerically, but is included as a guide 
to the eye since the chirality must vanish at /3 = 0. 



dictions of our theory (see Fig. HJi). At both transitions 
KA vanishes in a rather abrupt manner, which suggests 
that the transition is of the first order. One can also no- 
tice that the amplitude of the chiral order decreases as (3 
tends to 0, which is explained by the fact that ai /3 — 
the chirality should vanish since this limit corresponds to 
two decoupled unfrustrated ladders. 

For ferromagnetic rungs, J±/ J2 — —2, there is a tran- 
sition in Ks at positive /3 ~ 1.6 which has a similar be- 
havior to the corresponding transition in ka at > 0, 



while the situation at negative /3 is different: Looking 
simply at the chiral correlation functions, one tends to 
think that the chirality persists all the way up to /3 = —4, 
and merely the influence of the boundaries seems to in- 
crease considerably at /? < —3.0. Fig. [5] shows how the 
distribution of the magnetization along the tube changes 
in this region of (3. One can see that there is a transition 
at around /? ~ —3.2 which corresponds to the formation 
of a macroscopic magnon "droplet" in the ground state, 
sitting in the middle of the system. This is exactly the be- 
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FIG. 8. (Color online). The distribution of the rung magne- 
tization Mn = {S^^i + 5^,2) along the 5 = 1 spin tube of the 
length L = 128, at fixed J±/J2 ~ —2 and several different 
values of the frustration parameter /3 = J1/J2 in the vicinity 
of the transition from the repulsive to the attractive magnon 
gas. 
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FIG. 9. (Color online). The same as in Fig. |8] in the 
vicinity of the transition from the vector chiral into the two- 
component Tomonaga-Luttinger liquid phase. 



havior found in Ref . |l_9i for ferromagnetic frustrated S — 1 
chains, and indicates that the region with /3 < —3.2 ex- 
hibits a metamagnetic jump in the magnetization curve 
(the states with a "droplet" are never reahzed as true 
ground states at fixed magnetic field, they are only pos- 
sible if the number of magnons is artificially fixed) . 

It is worthwhile to remark that the behavior of the 
magnetization distribution can be also used to detect 
the transition between the VC and TLL2 phases. Fig. [S] 
shows how the magnetization oscillations, which are lo- 
calized at the boundaries in the VC phase, penetrate the 
bulk and spread over the entire system when one moves 
across the point /3 ~ 1.5. Comparing this behavior with 
Fig. El;, one can see that /3 « 1.5 is indeed the transition 
point where the symmetric chirality vanishes. 

For S — ^ tube, we have done similar calculations as 
for the spin-1 case, along two cuts at J± = ±2 J2 in the 
phase diagram. The resulting behavior of the chiral order 




frustration parameter p 

FIG. 10. (Color online). Chirality order parameters of the 
S — ^ tube extracted from the large-distance behavior of the 
correlation functions: (a) shows results obtained by means of 
the 5'f/(2)-symmetric DMRG method for L = 64 system (128 
spins) with J±/J2 = 2.0, in the sector with the total spin 
Stot = 116; the point at /3 = is included eis a guide to the 
eye; (b) the results obtained by the standard ?7(l)-symmetric 
DMRG calculation for L = 128 tube with Stot = 230 and 
J±/J2 = -2.0. 

parameters along those lines is shown in Fig.[TUl While at 
Jj_ > the picture is essentially similar to that for spin-1 
system, as described above, at negative J± the transition 
at lower values of /? looks rather different: the chiral order 
parameter Kg disappears in a very smooth way, as seen 
from Fig.llOb (the corresponding correlation functions of 
the symmetric chirality are presented in Fig. Illb ). At 
the same time, the distribution of magnetization at this 
transition shows the development of a spin density wave 
as seen in Fig. Illb . This is reminiscent of what happens 
in ferromagnetic frustrated spin chains,— i^^i^ and indi- 
cates that this transition corresponds to the formation 
of bound states of finite number of magnons, in contrast 
to the 5 = 1 case where there is a single bound state 
absorbing all the magnons present in the system. 

Comparing our numerical results for the selected cuts 
in the phase diagram with the analytical predictions of 
the two-component Bosc gas approach, one can see that 
our theory captures fairly well the physics of phase tran- 
sitions in the frustrated spin tube model ([1]). Compar- 
ing the results for the S = ^ spin tube (double zigzag 
chain) with those obtained by the same approach for sin- 
gle S = ^ zigzag chaiujiS one can see that the accuracy 
of the prediction for the VC-TLL2 transition is improved 
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FIG. 11. (Color online), (a) The distribution of the rung 
magnetization Mn = (5^,1 + S^,2) along the L = 128 spin-^ 
tube, at fixed J±/J2 = —2 and several values of the frustra- 
tion parameter P in the vicinity of the transition from the 
vector chiral phase to the phase with bound magnons and 
dominating spin density wave correlations. 

when one includes sufficiently strong rung coupling | J_l | . 
One obvious drawback of the theory, as mentioned at 
the end of Sect. IIIIl is that it does not take into account 
the possibility to have the lowest energy of the bound 
state at the total momentum different from ±2Q, which 
is realized for S = ^. We see that for that reason, our 
analytical predictions underestimate the size of the re- 
gion dominated by bound states for S = Apart from 
that, one may call the agreement between the analytical 
theory and numerical simulations satisfactory. 



we hope that our findings, which establish the high-field 
slice of the phase diagram, will stimulate experimental 
studies of field-induced phases in those systems. 

In the vicinity of saturation, the system can be repre- 
sented as a dilute gas of two flavors of bosonic quasipar- 
ticles corresponding to magnons with momenta around 
two degenerate incommensurate minima of the magnon 
dispersion. Using the method previously proposed for 
frustrated spin chains ^SiiiS we calculate effective inter- 
actions in this two-component Bose gas, and establish 
the high-field phase diagram of the frustrated spin tube. 
We show that the phase diagram contains two types of 
vector chiral phases (with symmetric and antisymmetric 
long-range chiral order), the two-component Luttinger 
liquid, and phases dominated by the presence of bound 
magnons. 

We complement our analytical results by the numerical 
studies of S' = 1 and S = frustrated tubes by means 
of the density matrix renormalization group technique. 
We analyze the behavior of chiral correlation functions 
and distribution of the magnetization along several cuts 
in the phase diagram, and extract the position of the 
corresponding phase boundaries. The numerical results 
are found to be consistent with our analytical predictions. 
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V. SUMMARY 

We study the ground state phase diagram of strongly 
frustrated four-leg spin-S* tube (which may be alterna- 
tively represented as two coupled zigzag spin chains, 
or as a two-leg spin ladder with next-nearest-neighbor 
couplings along the legs, see Fig. [T]) in a strong mag- 
netic field in the vicinity of saturation. The model is 
motivated by the physics of such frustrated quasi-one- 
dimensional spin-i materials as Sulfolane-Cu2Cl422ri2^ 
and BiCu2P06r^>2^"— but is also interesting in itself 
as the simplest model of coupled frustrated chains. Al- 
though both in SUI-CU2CI4 and BiCu2P06 the saturation 
field is too high to be accessible in current experiments. 



Appendix: Computing the Gross-Pitaevsky 
couplings Til, 

Consider in some detail the procedure of solving the 
Bethe-Salpeter equations (IT51) . To reduce the number of 
Fourier harmonics, we first symmetrize the kernel. In 
doing so, we use the identities = e"j,, V^^{k,k') = 

F"f(fc',fc), and rf{k,k';E) = T°fq{k' ,k] E). Let us 
introduce the following functions that are even in the 
transferred momentum q: 

Af ^Tf{Q,Q-^E,), (A.l) 
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then one has Tu = = 0) and Ti2{E) = B°'°'{q = 

Q), with a — a for J± > and a = c for J± < 0. One 
can rewrite Eqs. ([T8| as 



P 7 



PI ^ 



(A.2) 



where symmetrized kernels Ug^{p), v^^{p) are even func- 
tions of p: 



Assume for definiteness that Jj^ > 0, then of eight 
equations (jA.2p we need only a pair of equations for A^^ 
and , and another pair of equations for B'^^ and B^°' . 
Solutions to those equations can be now sought in the 



form containing only even Fourier harmonics: 

= X{) + xi cos q + X2 cos 2g, 
=yn+yi COS q + y2 cos 2g, 
+ xi cos g + 2^2 cos 2(j, 
yo + yi cos (7 + ^2 cos 2q. 



Bg 



(A.4) 



This ansatz transforms each of the above pairs of the 
integral equations into a system of 6 linear equations for 
6 variables. However, it follows from those equations that 



yi^xi, yi=xi, y2 = x2, y2 = x2, 



(A.5) 



so the size of the corresponding linear problems reduces 
to 4 X 4. For S = ^, one has to perform the limit U — 
+00. After solving the linear systems, one can read off 
the effective couplings 



Tii=xo+xi+X2, = a;o^ 



/3~ ^fl3' 



1)^2. (A.i 



The procedure for Jj_ < follows Eqs. (IA.4I) - (|A.6| ). with 
the obvious interchange a O c of the magnon branch 
labels. Below we list the equations for Xi, yi, Xi, yi in 
the form that is valid for any sign of J± as well as for 
any value of S, including S = ^. The resulting systems 
of equations can be cast into the following form: 
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where the matrix coefficients are given by 

1 r fd, 



I'- = 



^ Jo '^S{Jq+p + Jq-p — 2Jq) + 



1 



fifj 



fi = 1, /2 = cosp, fs = cos2p. 



r- = 



2S{Jq+p 

fifj 



Jq-p + 2\J^\ - 2Jq) + E, 



TT Jo 4S{\J±\ + Jp-Jq)+E,' 



(A.9) 



and can all be computed analytically in a closed, though 
somewhat cumbersome, form. The solutions for Fn, F12, 
which follow from Eqs. (IA.7I) . (lA.Sp . can be obtained with 
the help of any good computer algebra system (we used 



Maple), but are too bulky to be presented here (the re- 
sult, saved in a plain ASCII format, is several megabytes 
large) . 
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